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SUMM8RY 


fl  thecreticaS  study  of  the  fea’.Wy  of  the  moving  thermoacoustic 
source  es  a  sound  projector  in  underwater  acoustic  applications  is  presented. 

The  thermoacoustic  process  employs  direct  heating  of  the  acoustic 
medium  to  produce  a  controlled  local  thermal  expansion  which,  in  turn, 
generates  sound  wattes.  The  thermal  energy  is  delivered  by  a  laser  beam  to 
the  water  surface  without  any  physical  contact.  The  analysis  in  this  report  is 
limited  to  the  thermoacoustic  conuersion  process,  within  the  broader  class  of 
optoacoustlc  processes.  The  thermoacoustic  source  may  be  moved  by 
deflecting  the  laser  beam.  R  rotating  mirror  is  commonly  used  to  deflect  the 
beam.  The  deflection  causes  the  position  of  the  thermal  sound  source  to 
change  as  a  function  of  time.  This  allows  the  signal  to  be  Doppler  shifted  to  a 
degree  far  greater  than  that  achievable  with  conventional  sound  projectors 
that  hove  to  be  in  physical  contact  with  the  acoustic  medium. 

R  summary  of  the  findings  presented  in  this  report  1$  as  follows.  Since 
the  Ooppier  shift  is  direction  dependent,  it  was  found  that  the  moving 
optoacoustic  sources  are  potentially  useful  for  direction  fouling  by  virtue  of  tlie 
extremely  high  Ooppier  shifts  achievable.  While  the  thermoacoustic  array  has 
certain  advantages  over  conventual  acoustic  projectors,  notably  the 
noncontact  property  and  the  Doppler  direction  finding  capability,  its 
optoacoustic  energy  conversion  efficiency  is  no  better  than  that  of  the 
stationary  thermoacoustic  array.  This  conclusion  was  arrived  at  through 
mathematical  analysis  from  basic  principles,  and  supported  by  computer 


generated  numerical  examples.  The  conuersion  efficiency  mas  found  to  be 
strongly  dependent  on  the  acoustic  signal  carrier  frequency  and  on  the  optical 
signal  maueform.  Studies  of  the  effect  of  optical  input  signal  characteristics  on 
the  conuersion  efficiency,  the  efficiency  upper  bound,  and  the  pro jected  range 
penetration  under  oceanic  conditions  mere  made.  It  mas  found  that  the 
thermnacoustic  conuersion  mas  most  efficient  rnnen  the  optical  energy  mas 
defines  ed  as  an  impulse  train.  It  mas  concluded  that  the  thermoacoustic  array, 
mithin  the  broader  category  of  optoacoustic  processes,  is  expected  to  be  useful 
os  a  laboratory  tool  mhere  relatiuely  high  ultrasonic  frequencies  are  used,  but 
not  for  sonar  applications  because  of  the  lorn  thermoacoustic  conuersion 
efficiency  in  the  useful  range  of  sonar  frequencies.  The  efficiency  is 
fundamentally  limited  by  the  physical  properties  of  the  medium,  particularly  the 
coefficient  of  nermal  expansion  and  the  specific  heat  B  different 
optoacoustic  conuersion  process,  mhich  employs  nonlinear  physical  reactions  to 
generate  sound,  is  being  pursued.  It  promises  to  be  more  efficient  by  seuerai 
orders  of  magnitude. 


L  WTRODUCTION 


fl  “thermoacoustic"  source  Is  realized  when  heat  Is  directed  into  a  emafl 
region  of  an  acoustic  medium.  The  thermal  expansion  or  contraction  uf  the 
medium  causes  a  uolumetric  change  which  propagates  as  an  acoustic  mat* 
outwards  from  the  neated  region.  Since  the  thermal  source  is  often  generated 
optically,  it  is  sometimes  referred  to  as  a  “thermo-opticaP  source.  The 
optically  generated  thermoacoustic  source  also  belongs  to  the  broader  class  of 
“optoacoustic"  sources.  These  terms  may  be  used  interchangeably,  except  in 
those  cases  where  their  differences  are  important. 

bi  1976  the  process  was  investigated  theoretically  and  experimentally  by 
Muir,  Culbertson,  and  Clynch1  for  the  case  of  a  stationary  soiree.  The  acoustic 
energy  produced  by  the  thermoacoustic  array,  (Tfl),  with  any  reasonable 
quantity  of  optical  energy  through  the  thermoacoustic  conversion  process,  was 
found  to  be  small  It  was  several  orders  of  magnitude  smaSer  than  that 
produced  by  conventional  acoustic  soirees  under  equivalent  operating 
conditions. 

It  is  stated  by  Lyamsheu  and  Naugol'nykh^  in  a  comprehensive  review 
paper  that 

Under  the  conditions  of  thermo-optical  generation  of  sound  (Le,  the 
thermal  mechanism)  the  conversion  efficiency  attains  at  best  10*^  -  10~6." 

Therefore,  f or  underwater  detection  ond  communication  purposes,  the 

thermoacoustic  conversion  efficiency  must  be  significantly  increased  before 

the  TR  can  be  considered  as  a  viable  acoustic  source.  There  have  been  several 

recent  Soviet  papers  on  the  subject  with  diverse  claims,  particularly  with 

regard  to  the  moving  thermoacoustic  source. 
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It  was  suggested  by  some  that  the  moulng  thermoacoustJc  source,  (MTfl), 
can  giue  significantly  higher  conuersion  efficiencies.  The  simptest  case  is  that 
of  an  MTfl  traueling  in  a  straight  line.  Here  Bozhkou  et  al,3  claim  that  mouing 
thermo-optical  sources  are 

highly  directional  sound  sources  with  a  tunable  frequency  and  considerably 
higher  efficiency  than  that  of  stationary  sources." 

in  an  earlier  publication,  Bozhkou  et  at4  also  claimed  that 

"the  3,uen  method  of  sound  generation  affords  a  realistic  possibility  for  the 
creation  of  powerful  short  pulses  with  a  high  frequency  carrier." 

r  e  case  of  an  MTfl  executing  uariou*  types  of  motion,  including 
oscillatory  motion,  has  been  studied  by  Lyamsheu  and  Sedou,3  but  they  clair., 
only  that  the  efficiency  of  mouing  sources  is 

practically  the  same  as  for  a  stationary  source." 

It  is  agatr<  t  this  background  of  diuerse  and  conflicting  claims  regarding 
the  perfo?  jnce  of  the  MTfl  that  the  inuestigatlon  takes  place. 

Rn  oueruiew  of  the  report  is  briefly  as  follows. 

The  thermoacoustic  conuersion  process  is  analyzed  from  basic  principles 
in  section  3.  Both  the  stationary  and  mouing  cases  are  examine-i  The  special 
case  of  a  Mach  number  of  unity  is  also  treated.  The  conuersion  efficiencies 
ore  compared  in  section  4.  fhe  influence  of  the  optico1  signal  wnueform  and 
the  acoustic  signal  frequency  on  the  conuersion  efficiency  is  examined  and  then 
the  energy  conuersion  efficiency  upperbountis  are  established  In  section  5.  The 


performance  of  the  thermoacoustically  generated  signal,  particularly  range 
penetration  under  typical  oceanic  conditions,  is  analyzed,  supported  by 
numerical  examples,  in  section  6.  The  direction  finding  capability  is  discussed  in 
section  7.  The  findings  are  discussed  and  compared  with  the  claims  quoted 
aboue  in  sections  8  and  9.  Finally,  the  remaining  critical  issues  and  the  methods 
to  resolue  them  are  discussed  in  section  10. 


» 


!L  OSJECTIUf ; 

1.  The  thermoacoustic  process  and  its  conuersion  efficiency  Luill  be 
investigated  from  first  principles. 

2.  The  upper  limit  of  the  conuersion  efficiency  and  the  factors  which 
affect  it  will  be  determined 

3.  The  claims  that  have  been  made  for  the  thermoacoustic  process  in  the 
open  literature  will  be  examined 

4.  The  inability  of  the  thermoacoustic  source  as  a  projector  of  acoustic 
signals  m  underwater  acoustic  applications  will  be  assessed 
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Hi.  THEORY 


R.  Basic  Theory 

In  this  study,  the  medium  is  water  and  the  means  of  introducing  heat  is  a 
laser.  The  thermal  input  may  be  controlled  by  modulating  the  laser  power.  The 
absorption  coefficient  of  light  in  water  is  assumed  to  be  uniform.  Therefore 
the  thermoacoustic  source  is  exponentially  tapered  as  a  function  of  penetration 
distance  into  the  water.  The  geometry  is  shown  in  Fig.  1. 

R  suitable  starting  point  is  the  expression  for  the  Fourier  transform  of 
the  far  field  acoustic  pressure  generated  by  a  stationary  TR  such  as  the  one 
derived  by  Lyamshev.^  Lyamsheu's  expression  was  for  a  lossless  medium  but 
in  the  following  expression,  which  is  derived  in  Appendix  ft,  the  absorption 
coefficient  is  included. 

PaoM  =  k  CQ  Dq  expt  ikr0  -  ar0 )  /r0 ,  (1) 

where 

Ca  -  RHc/HTTCp)  (2) 

and 

k  =  w/c .  (3) 

The  parameters  Ca  and  Da  haue  the  following  significances.  Ca  *  presents  the 
thermoocoostic  conversion  factor  afforded  by  the  physical  properties  of  the 
medium.  0Q  represents  the  acoustic  losses  due  to  diffraction  taused  by  the 


9 


( 


AIR 


ftaite  array  length  and  width,  and  the  optical  loss  due  to  reflection  at  the 
surface.  In  Edition  to  the  general  case,  the  special  case  inhere  the  optical 
beam  is  uertical  and  the  intensity  profile  in  the  horizontal  plane  is  a  Gaussian 
function  is  worth  considering  because  it  is  capable  of  yielding  a  dosed  form 

solution;  in  this  case,  Da  is  gioen  by2 

Da  =  jikcosR  exp{  -k2a2sin28  /4  V(u2  +  k2cos2(6) ) ,  (4) 

where  a  is  the  1/e  radius  of  the  6aussian  functioa 

The  response  to  an  optical  impulse  is  easily  obtained  by  equating  U)Q((0) 
to  the  Fourier  transform  of  an  impulse. 

B.  The  Mouinq  Thermoacoustic  flrrau 

Let  us  now  consider  the  case  of  on  MTR.  The  MTfl  wiU  be  mrdeied  as  a 
continuum  of  stationary  optical  impulses  distributed  along  its  path  in  space  and 
time.  Let  the  water  surface  be  the  x-y  plane  in  a  Cartesian  coordinate  system, 
as  shown  in  Fig.  1.  The  top  end  of  the  MIR  must  always  be  at  the  water 
surface.  Let  the  changes  in  its  x  and  y  coordinates  as  a  function  of  time  be 
described  by  path  functions  x(t)  and  y(ti,  where  the  time  t  is  the  independent 
uariable,  and  let  the  path  be  centered  at  the  origin.  Therefore,  the  response  of 
an  MTR  of  arbitrary  path  function  and  an  optical  power  input  function  w(t)  may 
be  formulated  as  a  connotation  of  the  impulse  response  and  the  source 
distribution  function  in  space  and  time.  Since  time  is  the  independent  uariatrie, 


11 


the  result  may  be  expressed  as  the  time  integral 


Pa(<*>)  *  k  Cg  .^J00  Da  uMt)  e»p(ik[r  -  ctj  -ar)  }x  dt,  (5) 

where  r,  as  shown  in  Fig.  1,  is  the  distance  from  the  surface-end  position  of 
the  MTR  oHx(t),y(t),0)  to  the  receiuer  position  at  (x^O.z-j).  It  is  implicit  that  r 

uaries  as  a  function  of  t  The  acoustic  pressure  pa  may  be  obtained  by  taking 
the  inuerse  transform 

Pa(t)  =  _mr  P3«o)  exp(-jwt)  dw .  (6) 

Using  Parseual's  theorem,  the  acoustic  energy  density  Ea  is  giuen  by 

Ea  =  -ooJ“iPaMi2/(2hpc)dW.  (7) 

C.  The  Stationa> .  "termoacoustic  Array 

In  the  stationary  case,  the  source-receiuer  distance  rQ  is  a  constant. 

This  allows  Eq.  (5)  to  be  considerably  simplified,  gluing 

PaM  *  ( kCaDg  /r0)  expiikt  0-ar0)  _M|°®w0(ttexp(-ikct)dt .  (8) 

The  integral  in  Eq.  (8)  is  recognized  os  a  Fourier  integral  Therefore  the  result 
may  be  expressed  in  terms  of  the  Fourier  transform  of  the  optical  power 
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function  1D0M 

PaUo)  =  (kCaOa/r0)  exp(ikr0 -arQ)  LU0M  5  (9) 

which  is  identical  to  Eq.  (1).  This  confirms  that  the  expression  for  the  acoustic 
signal  m  Eq.  (5)  is  uatid  for  both  mouing  and  stationary  cases. 


* 


ttJ.  THE  EFFECT  OF  MOUEMENT  ON  THE  OPTICAL  ENERGY  REQUIREMENT 

To  examine  the  effect  of  mouement  on  the  optical  signal  requirements, 
let  us  turn  Eq.  (5)  around  and  keep  the  acoustic  signal  constant.  Assuming  that 
the  signal  may  be  uniquely  defined  by  its  Fourier  transform,  this  problem  is 
equiualent  to  that  of  examining  the  factors  affecting  the  Fourier  transform  of 
a  giuen  acoustic  signal.  The  energy  requirement  analysis  mill  start  Luith  the 
long  ranye  case  where  the  MTA  is  mooing  at  a  constant  Mach  number  relatiue 
to  the  receioer.  The  term  'long  range"  is  used  to  indicate  the  situation  where 
source-receioer  distance  is  such  that  the  launch  and  arriual  angles  of  the  sound 
rays  do  not  change  significantly  ouer  the  duration  of  the  signal.  Then,  the 
results  will  be  generalized  to  include  the  case  of  uarying  Mach  number.  Finally, 
the  long  range  requirement  is  relaxed.  In  all  cases,  the  receioer  is  assumed  to 
be  in  the  farfield  of  the  source. 

Let  us  consider  the  long  range  case  where  the  MTR  is  traoeling  at  a 
constant  Mach  number  relatiue  to  the  receiuer.  Referring  to  Fig.  1,  the  mean 
position  of  the  MTR  is  centered  at  the  origin  qf  a  Cartesian  coordinate  system, 
and  the  receiuer  is  located  in  the  x-2  plane.  It  is  assumed  in  the  long  range 
case  that  the  MTA  path  length  is  small  compared  to  the  source-receioer 
distance;  therefore  changes  in  the  source-receioer  distance  and  the  depression 

angle  are  negligibly  small  compared  to  their  mean  uaiues,  r0  and  8^  It  also 

follows  that  the  y  component  of  the  MTA  ueiocity  has  a  negligible  effect  on  the 
Mach  number.  In  this  case,  Eq.  (5)  may  be  approximated  by 
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(10) 


Pa(<o)  =  (kCa0a/ro)  exp(ikr0 -ar0) 

-mr  iu(t)  exp(ik(uHt  cos(80)  -  ct] )  dt, 

where  ux  is  the  component  of  the  MTfl  velocity  in  the  x  direction.  The  integral 

in  Eq.  (10)  is  recognized  as  a  Fourier  integral  and  it  mag  be  expressed  as  the 
Fourier  transform  of  w(t),  gluing 

Pq«a))  =  (kC0D0/r0)  ero(ikr0  -ar0)  LU(a)(1-M)) ,  (11) 


where 

M  =  (ux/c)cos  (0) . 

fl.  The  MM1  Case 

Comparing  Eqs.  (9)  and  (11)  for  the  case  luhere  M#1,  the  differe  xe 
between  the  optical  power  signals  in  the  stationary  and  moutng  cases  can  be 
cleerlg  seen.  In  order  to  obtain  the  same  acoustic  signal,  the  fatter  must  be  a 
Doppler  shifted  uersion  of  the  former, 

I1)(W(1-M))  *  Uf0M.  (12) 

If  the  total  optical  energy  expended  to  produce  a  giuen  acoustic  signal 
can  be  somehow  reduced  through  mouement  of  the  thermocaoustic  array  then 
there  would  be  a  gain  in  conuersion  efficiency.  Consider  the  totol  optical  energy 


which  is  giuen  by 


WQ)  =  dt.  (13) 

The  total  optical  energy  is  recognized  as  the  Fourier  transform  of  the  optical 
power  at  (*>=0.  Since  a  Doppler  shift  influences  the  Fourier  transform  at  all 
ualues  of  (*)  except  at  00=0,  there  should  be  no  difference  in  the  Fourier 
transforms  at  co-0,  that  is, 

ID(0)  =  UJo(0).  (14) 

Substituting  from  Eq.  (14)  into  (13),  it  is  clear  that  motion  has  no  effect  on  the 
total  optical  energy  requirement, 

w(t)  dt  *  J\(t)  dt.  (15) 

The  result  may  be  extended  to  the  uarying  Mach  number  case,  where 
the  x-coordina*e  of  the  MTfl  is  generally  represented  by  the  function  x(t).  The 
Fourier  transform  of  the  acoustic  signal  Is  now  "tuen  by 

Pa(oo)  -  ( k  Ca  Da  /r0 )  exp(ikr0  -ar0) 

^oJ^wd)  exp(itfxii)  cos(8p)  -  ctl )  dt.  (16) 

In  this  case,  the  optical  power  function  w(t)  needs  to  be  nonuniformly  Doppler 
shifted  to  maintain  the  same  acoustic  signal  If  the  nonuniform  Doppler  shift  is 


1? 


one  that  changes  smoothly  as  a  function  of  time,  then  it  may  be  approximately 
broken  down  Into  a  number  of  discrete  time  interuals  in  which  the  Doppler 
shifts  are  uniform, 

Pa(«)  «  (kCaDa/r0)  exp(ikr0 -ar0) 

N 

.oJ00  2wn(t)  exp(ik|uxnt  cos(80)  -  ctl  j  dt,  (17) 
n=1 

where 

iun(t)  =  w(t)  if  tn<  t  <tn+1 
0  otherwise 


and 

uxn  =  ^n+1* '  x(tn))/(tn+i  -  tn) . 

The  approximation  becomes  more  exact  as  the  number  of  time  interuals  tends 
to  infinity.  Since  it  has  been  concluded  that  uniform  floppier  shifts  do  not 
affect  the  optical  energy  requirement,  it  follows  that  the  total  optical  energy 
requirement  summed  ouer  all  the  frne  interuals,  where  the  floppier  shift  in 
each  interuaJ  is  uniform,  must  also  be  unaffected. 

N 

x  uyo)  -  uyo) ,  <i8) 

n*1 

Therefore,  Eq.  (15)  remains  uaiid  and  the  optical  energy  requirementremains 
unchanged. 


Finally,  the  analysis  may  be  broadened  to  include  the  general  farfieid 
case  where  changes  in  r  and  8  cannot  be  considered  as  small  Equation  (5)  may 
be  rewritten  in  a  form  similar  to  Eq.  (16), 

Pa«*»  «  ( k  Ca  Dac  /r3 )  exp(ikr0 -ar0) 

w(t)  (Da/Dao)  exp(ik[r  -  -  ctl  -a(r  -  r0))  dt .  (19) 

In  this  case,  the  optical  power  function  wQ(t)  must  not  only  be  distorted  to 
compensate  for  the  changing  Doppler  shift,  it  must  also  be  further  distorted  to 
compensate  for  the  differential  changes  in  the  diffraction  loss,  (Da/Dao),  and 

absorption  loss,  exp(-a(r  -  rQ)).  Since  the  losses  are  range  dependent  and  the 

range  may  change  significantly  ouer  the  path  of  the  MTR,  the  comparison 
between  the  mouing  and  stationary  cases  is  not  straightforward.  Howeuer, 
since  r0  and  80  are  mean  ualues  of  r  and  8  ouer  the  significant  part  of  the 

path,  the  net  increase  or  reduction  in  the  required  optical  input  to  maintain  the 
same  acoustic  signal  must  be  a  second  order  quantity.  Therefore,  the  optical 
energy  requirement  is  unchanged  at  least  to  the  first  order  of  accuracy. 

B.  The  M-1  Case 

let  us  now  examine  the  behauior  of  the  MTR  when  the  Mach  number  M 
is  equal  to  unity.  From  Eq.  (11),  it  can  be  seen  that  the  ualue  of  the  IDM1-M)), 
for  all  finite  ualues  of  Q,  is  equal  to  LU(O).  That  is,  the  optical  excitation 


... 1  wy  '  mm  >-mj ,  <— y 


appears  to  collapse  into  an  impulse.  Consequently,  in  the  absence  of  any 
distortions,  the  acoustic  maue.orm  mould  be  the  deriuatiue  of  an  impulse.  The 
acoustic  signal  is  obtained  by  setting  M=1  in  Eq.  Ml), 

P8(%1  -  ( k  Ca  Da/rQ  )  exp(ikr0 -ccr0)  UMO) .  (20) 

Note  that  the  signal  is  dependent  only  cm  the  optical  energy  IlHO). 
Therefore,  the  shape  of  the  optic?!  pomer  maueform  has  no  influence  on  the 
acoustic  output  The  acoustic  maueform  is  only  determined  by  the  frequency 
dependence  of  the  diffraction  loss  Da,  mhich  is  controlled  by  the  physical 

dimensions  of  the  thermoacoustic  source,  and  the  absorption  coefficient  a. 
Homeuer,  this  does  not  mean  that  the  thermoacoustic  conuersion  efficiency  of 
an  MTfl  traueiing  at  Mach  1  mill  be  any  higher  than  that  of  a  stationary  TR. 

In  the  case  of  a  stationary  TR,  an  identical  result  can  be  obtained  by 
using  an  optical  impulse.  In  this  cose,  the  Fourier  transform  of  the  optical 
pomer  U)q(cj)  mill  be  equal  to  U)0(Q)  for  afl  significant  ualues  of  (*).  Substituting 
into  Eq.  (9), 

^V^impulse  *  ^  ^a  ®a^ro  ^  exp(ikrQ  -ajro)  IDq(O)  .  (21) 

From  Eqs.  (20)  and  (21),  it  is  euident  that  if  the  optical  energies  are  equal, 
i.e.,  ll)(0)*  iDo(0),  then  the  Fourier  transforms  of  the  acoustic  signals,  hence  the 

acoustic  moueforms  themselues,  must  be  identical.  Consequently,  it  is 
concluded  that  a  stationary  TR  energized  by  an  optical  impulse  mould  produce 


exactly  the  same  acoustic  signal  as  an  MTfl ,  mouing  at  M=t  n  the  direction  of 
the  receiuer,  for  the  same  optical  energy.  The  signals  mill  be  identical  in  both 
waueform  and  energy  density. 

Substituting  from  Eq.  (20)  or  (21)  into  Eq.  (7),  an  expression  for  the 
acoustic  signal  energy  (tensity  may  be  obtained, 

Ea  =  C02  Girpcf1  r0'2  llXO)2  k2  0o2  exp^or^  dm .  (22) 

Note  that  the  Fourier  transform  of  the  optical  power  input  may  be 
brought  outside  the  integral  because  it  is  constant  for  aB  significant  ualues  of 
frequency.  The  expression  also  shows  that  if  the  absorption  coefficient  a 

mere  finite  and  independent  of  frequency,  and  the  Effraction  loss  Da  increases 
with  frequency  by  a  rate  less  than  the  second  power,  then  the  acoustic  energy 
density  Ea  mould  be  inf inite.  to  practice,  however,  a  monotonicafly  increases 
with  frequency  and  this  alone  is  sufficient  to  ensure  that  the  acoustic  energy 
density  is  always  bounded,  to  practice,  the  diffraction  loss  DQ  also  increases 

rapidly  with  frequency  beyond  a  cutoff  frequency  determined  by  the  physical 
dimensions  of  the  array.  There  are  also  other  losses  not  represented,  such  as 
that  due  to  finite  amplitude  effects,  which  would  further  limit  the  acoustic 
energy  density  at  the  recemer. 
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R,  Upper  Bound  of  the  Energy  Density 

Since  it  has  been  demonstrated  in  the  foregoing  sections  that 
thermoacoustic  conversion  by  stationary  and  mouing  sources  are  equiualent, 
the  upper  bound  def  ined  below  will  be  applicable  to  both  cases.  It  is  possible  to 
obtain  &n  upper  bound  for  the  ualue  of  the  energy  density,  Ea,  by  using  the  fact 

that  the  term  0a  in  Eq.  (51  and  subsequent  equations  is  a  loss  term  and  its  ualue 
is  always  less  than  unity.  Thus,  setting  Da  equal  to  unity  in  Eq.  (9)  and 

substituting  into  Eq.  (71,  the  upper  bound  may  be  expressed  by  the  following 
inequality, 

Ea  <  (Ca/r0)2(2npcf1  k2e>:p(-2  ctr0J  |lil(u)(1-MI)j2  dw  .  (23) 

It  can  be  seen  that  the  Fourier  transform  of  the  optical  signal, 

UI(cj(1-M)I,  at  nonzero  ualues  of  the  acoustic  signal  frequency,  (  w*0  ),  has  a 
direct  influence  on  the  acoustic  energy  density  but  no  direct  connection  with 
the  optical  energy  input.  On  the  other  hand,  the  ualue  of  UHwtl-M))  at  u=0, 
which  represents  the  total  optical  energy,  has  no  direct  effect  oh  the 
generation  of  the  acoustic  signal.  Since  these  two  Fourier  transform 
components  are  mutually  exilusiue,  it  is  conteiuable  that  an  acoustic  signal  can 
be  generated  without  any  net  optical  energy  dissipation.  Howeuer,  in  practice, 
this  is  not  achievable  because  negutiue  thermal  power  dissipation,  i.e.,  the 
absorption  of  thermal  power  from  the  water  by  the  laser  beam,  is  not  possible. 
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The  uoJue  of  llKoMl-M))  at  0)^0  represents  a  positiue  energy  teas  which 
attorns  the  t  her  mat  power  changes  necessary  for  acoustic  generation  to  take 
place  without  reuersing  the  direction  of  optical  pouter  flout  in  the  laser  beam. 

Under  this  restriction,  it  can  be  shoutn  that  the  ualue  of  UKq(1-M))  at 
u*0  cannot  be  greater  than  the  uaiue  at  ot=0,  Le^ 

<  UHO)  for  (0*0.  (24) 

The  proof  is  giuen  in  Rppendix  B.  Substituting  from  En.  (24)  into  Eq.  (23),  the 
acoustic  energy  density  upper  bound  is  obtained  in  terms  of  the  optical  energy, 
11K0), 

Ea  <  (Ca/r0)2(2Tfpc)_1  UKO)2  .cj"  k2exp<-2  arj  dot .  (25) 

The  upper  bound  giuen  by  Eq.  (25)  is  uoitd  in  ail  cases,  including  the  M-1 
case.  For  narroutbond  signals,  the  upper  bound  may  be  more  exactly  expressed 
in  terms  of  the  mean  uiaue  number,  kQI  and  the  acoustic  bandwidth,  Bg  ( in  Hz), 

£a  narrcuband  <  «Vro»2'''Pcl»W0>2  K02exp<-2  ar„)  2Ba .  (26) 

For  a  giuen  bandwidth,  Ba,  and  range,  xv  there  is  an  optimum  frequency  at 

which  the  narrowband  upper  bound  is  at  a  maximum.  Using  the  absorption 
curue  for  seawater  at  10°C,  a  plot  of  the  optimum  frequency  as  a  function  of 
range  was  compute*  shown  in  Fig.  2.  The  signal  bandwidth  was  assumed  to 
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FIGURE  2 

FREQUENCIES  AT  WHICH  UPPER  BOUNDS  OF  ENERGY  DENSITY 
AND  SIGNAL-TO-NOISE  RATIO  ARE  MAXIMIZED 
THE  NOISE  IS  ASSUMED  TO  BE  KNUDSON  NOISE, 

WHICH  INCREASES  AT  17  dB  PER  DECADE  OF  FREQUENCY 
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be  a  constant 


B,  Waveform  and  the  Energy  Conversion  Efficiency 

In  Eq.  (23),  it  is  evident  that  the  acoustic  signal  energy  density  is 
generated  by  the  square  of  the  magnitude  of  the  optical  power  Fourier 
transform,  UJMI-M)).  Therefore,  the  conversion  efficiency  of  waveforms  is 
directly  governed  by  the  relative  value  of  UKwO-M))  within  the  signal  band  to 
the  optical  energy  liHO).  For  a  given  acoustic  signal,  the  conversion  efficiency 
Is  maximized  by  mamtainmg  the  desired  values  of  W(w(1-M))  while  minimizing 
the  optical  energy  bias  liftO)  as  much  as  possible.  The  most  efficient  optical 
signal  is  one  that  achieues  the  equality  in  Eq.  (24)  over  the  significant  part  of 
the  signal  band,  in  Table  1,  the  peak  values  of  Wfcod-M))  within  the  signal  band 
are  compared  for  three  representative  waveforms;  in  each  case  UKO)  is  set 
equal  to  unity.  The  waveforms  are  the  sine  wave,  the  square  wave,  and  the 
impulse  train  as  shown  in  Fig  3.  In  each  case,  it  is  assumed  that  the  signal  band 
is  centered  on  the  fundamental  frequency  tbg. 

It  is  clear  that,  of  the  three  waveforms,  the  impulse  train  produces 
the  maximum  acoustic  signal  energy  for  the  same  optical  energy.  It  appears  to 
oe  the  most  efficient  waveform  for  thermoacoustic  generation  in  geiieral  as 
long  as  unwanted  harmonics  do  not  interfere  with  the  desired  signal. 
Furthermore,  it  is  well  suited  to  laser  technology  where  Q-switching  methods 
for  producing  high  energy  impulses  ore  weil  developed. 
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TR8LE I 

THE  FUNOflMENTRL  COMPONENTS  OF  CERTfflN  PERIODIC  UMJEFORMS 


Waveform 


MXw0(1-M)i!2  at  tlKO)  »1 


sine  wave  0-25 

square  wave  0.41 

impulse  train  1.00 


C.  Maximization  of  the  SiqnaHo-Noise  Ratio 


Consider  an  ideal  situation  where  there  are  no  boundaries,  the  speed  of 
sound  is  constant  and  isotropic,  the  background  noise  is  isotropic,  and  the 
receiver  employs  a  matched  filter  to  maximize  the  signaHo-noise  ratio  (S/N). 
For  a  one  way  transmission  between  on  acoustic  source  and  receiver,  the  S/N 
may  be  simply  expressed  by  the  fofloumg  equation . 

S/N  *  10tog(Ea)  -  lOlogtF^  +  .  (27) 

The  contribution  from  the  acoustic  source  in  this  equation  is  summed  up  in  the 
acoustic  energy  density  term,  Eg.  For  a  receiver  of  a  predetermined  directivity 
index  or  for  an  omnidirectional  receiver,  the  S/N  is  maximized  by  maximizing 
the  ratio  Eg/F^  For  narrowband  signals,  there  is  a  carrier  frequency  at  which 

the  upper  bound  S/N  is  a  maximum.  Using  the  expression  for  the  upper  bound 
signal  energy  density  m  Eq.  (26),  and  assuming  that  the  receiver  is  in  a  Knudson 


noise  limited  situation  ivhere  the  equivalent  plane  rnaue  noise  spectrum 
Increases  at  a  rate  of  1?  tiB  per  decade,  the  carrier  frequency  at  mhich  the  S/N 
upper bound  is  maximized  is  shoivn  as  a  function  of  range  in  Fig.  2. 

D.  Conversion  Efficiency  Uooer  Bound 

Let  us  define  efficiency,  n,  as  the  acoustic  energy,  Ul^  produced  per  unit 

of  optical  energy  input,  18(0),  Since  the  sound  is  assumed  to  be  ratiated 
dommvards  from  the  surface,  the  acoustic  energy  upper  bound  is  equal  to  the 
energy  density  giuen  by  inequality  (251,  without  the  absorption  loss  term, 
integrated  over  a  hemispherical  boundary  beKu  the  water  surface, 

I8a  <  (Cg2/(pcj)  UKO)2  _oJOT  k2  do .  (28) 

The  efficiency  upper  hound  is  then  found  by  dividing  the  acoustic  energy  by  the 
optical  energy, 

D  <  (Ca2/(pc))  18(0)  _mS°°  k2  d(o  .  (29) 

From  the  above  inequality,  it  appears  that,  in  the  absence  of  diffraction  loss 
and  any  other  type  of  losses  internal  to  the  thermoacoustic  source,  the 
efficiency  upper  bound  is  infinite.  This  is  clearly  erroneous  since  it  contradicts 
the  latv  of  conservation  of  energy.  The  flam  lies  in  the  assumption  that  the 
medium  reacts  to  the  application  of  heat  by  expanding  instantaneously.  LUhile 
this  assumption  is  valid  over  the  band  of  detectable  acoustic  frequencies,  it  uriB 


break  down  beyond  a  cutoff  frequency  determined  by  the  tag  time  of  the 
thermal  expansion  process.  The  cutoff  frequency  win  effectively  make  the 
limits  of  the  integral  in  P+  tz9J  finite  and  constrain  the  efficiency  upper  bound 
to  a  uaiue  less  than  or  equal  to  one. 

From  the  above  analysis,  it  can  be  seen  that  the  efficiency  upper  bound 
varies  over  a  wide  range  of  values,  depending  on  acoustic  signal  frequency.  It 
is  seen  that  for  high  frequency  sonar  te.g.,  30  kHzJ  and  reasonable  optical 
energy  levels  te.g.,  100  Jl,  the  efficiency  upper  bound  is  around  10  6,  which  is 
marginally  useful  The  range  penetration  of  a  signal  in  this  category  is 
discussed  later,  flt  long  range  sonar  frequencies  te.g.,  3  kHz),  and  at  reasonable 
levels  of  optical  energy  (e.g.,  100  J),  the  efficiency  upper  bound  is  less  than 
10'9  land  even  with  1000  J  of  optical  energy  it  is  still  less  than  10~8i 

For  narrowband  signals  the  efficiency  is  significantly  lower  than  that 
allowed  by  the  upper  bound,  since  efficiency  increases  with  bandwidth.  The 
expressions  for  total  acoustic  energy  and  efficiency  are  a  little  simpler  since  it 
Is  reasonable  to  approximate  the  value  of  the  wave  number  by  a  constant  over 
the  passband.  The  expression  for  the  total  acoustic  energy  upper  bound  is 
obtained  from  inequality  (281  by  replacing  the  integral  with  a  simple  product, 

narrowband  <  UJfQt^  k^  4irBa  (30) 

The  efficiency  upper  bound  is  then  given  by 

^narrowband  <  ^a^Pc^  ^  ^^a  •  ^1) 


As  the  kQ2  term  suggests,  the  gains  are  concentrated  at  the  high 

frequency  end  of  the  spectrum.  In  the  laboratory,  where  oper  ating 
frequencies  are  generally  higher  and  propagation  distances  shorter  than  in  the 
sea,  the  Tfl  promises  to  be  a  uery  efficient  and  useful  tool7  In  sonar 
applications,  howeuer.  absorption  hmU  uery  rapidly  eliminate  the  high 
frequency  components  before  any  useful  distance  can  be  couered.  The 
efficiency  upper  bounds  of  a  Tfl  at  certain  ualues  of  frequency  and  optica! 
energy,  assuming  that  the  bandwidth  Bg  is  one-half  of  the  carrier  frequency, 
are  shown  in  Fig.  4. 


Ul  NUMERICAL  EHIMT.ES  OF  MOtflNG  THERMORCOUSTIC  RBRflYS 


The  f  otowing  numerical  example  wifl  be  used  to  set  the  main 
characteristics  of  the  MTfl  into  (roper  perspective.  For  this  purpose,  an  MTfl 
ivas  modeled  for  the  following  set  of  operating  conditions : 

The  heat  source  urns  a  laser  directed  verticefly  downwards.  The  laser 
beam  profile  was  Gaussian  shaded  with  a  1/e  radius  of  8  mm.  The 
laser  defivered  a  total  of  190  J  of  optical  energy  at  a  wavelength  of 
1.06  jim,  in  a  train  of  ten  ptdses  at  a  repetition  rate  of  2  kHz.  The  MTfl 
traueled  at  M=0.9  along  the  x  axis  giuing  the  acoustic  signal  a  Doppler 
shifted  frequency  of  20  kHz  in  the  forward  (+x)  direction.  The  water 
was  assumed  to  be  seawater  at  a  temperature  of  10°C  with  an 
isotropic  noise  spectrum  equivalent  to  that  of  sea  state  4.  The  acoustic 
signal  was  received  by  an  omnidirectional  hydrophone  equipped  with  a 
matched  filter. 

The  energy  density  pf  the  signal  as  a  function  of  range  and  direction  was 
computed.  The  one-way  slant  range  penetration  contour  for  an  S/N  at  the 
receiver  of  20  dB  is  shown  on  Fig.  5.  The  computations  showed  that,  at  a 
depression  angle  of  8°  in  the  forward  direction,  the  range  penetration  reaches 
a  maximum  value  of  1400  m.  The  acoustic  source  level  was  computed  to  be 
165  dB  re  1  jiPa. 

By  numerically  integrating  the  energy  density  of  the  acoustic  signal  over 
a  hemispherical  boundary,  the  total  energy  conversion  efficiency  was  computed 
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to  be  1  x  IQ'8.  From  Eq.  (31),  the  upper  bound  of  efficiency  at  a  carrier 
frequency  of  zO  kHz  and  a  bandwidth  of  2  kHz  is  estimated  to  be  1.6  x  10'7 
Therefore,  the  actual  efficiency  is  considerably  less  than  the  upper  bound.  The 
discrepancy  is  due  to  the  reduction  in  carrier  frequency  as  a  function  of  angle 
away  from  the  forward  direction,  and  to  diff ration  loss. 

The  range  penetration  of  a  stationary  Tfl  pulsed  at  a  frequency  of 
20  kHz,  operating  under  the  same  conditions  and  of  equal  optical  energy,  is 
shown  in  Fig.  6.  It  is  seen  that  the  stationary  Tfl  has  the  same  acoustic  signal 
frequency  of  20  kHz  and  achieues  the  same  range  penetration  as  the  MTR.  In 
the  stationary  case,  the  acoustic  output  is  independent  of  azimuth.  It  has 
greater  coverage  but  no  directiuity.  The  total  energy  efficiency  was  computed 
to  be  1.3  x  IQ"7.  The  discrepancy  between  the  computed  efficiency  and  the 
upper  bound  of  1.6  x  iq"7  is  due  to  diffraction  loss.  Rlthough  the  stationary 
Tfl  achieues  better  oueraH  efficiency  than  the  MTfi,  the  gain  is  in  direction? 
other  than  the  forward  direction.  The  level  of  the  20  kHz  signal  in  the  forward 
direction  is  unaffected, 

Rn  indication  of  the  range  penetration  Gf  the  acoustic  signal  as  a 
function  of  its  frequency  can  also  be  extracted  from  Fig.  5.  The  direction  and 
range  penetration  of  the  10  kHz  and  the  5  kHz  signals  we  shown.  Since  the 
range  penetration  has  been  shown  to  be  the  same  whether  die  thermoacoustic 
source  is  moving  or  stationary,  the  following  conclusion  is  appficabie  fo*  «*.ih 
coses.  The  results  suggest  that,  for  the  given  optical  energy  of  190  J,  the 
range  penetration  below  20  kHz  decreases  as  the  frequency  is  reduced,  in  spite 
of  the  lower  absorption  toss  at  lower  frequencies.  This  i$  Demise  the 
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FIGURE  6 
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reduction  in  efficiency  (plus  the  increase  in  ambient  noise)  outweighs  the 
reductions  in  absorption  loss.  Further  numerical  analysis  showed  that  the 
range  penetration  peaked  at  20  kHz,  and  that  beyond  20  kHz  it  declined  with 
increasing  frequency.  Therefore,  it  is  concluded  that  at  an  energy  ieuet  of 
190  J  (or  less)  a  thermoacoustic  source,  mooing  or  stationary,  will  not  prouide 
the  range  penetration  required  for  medium  or  long  range  communications. 
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m.  DIRECTION  FINDING  CRPRBILITV 


In  spite  of  its  inefficiency,  the  MTH  has  a  potentially  useful  but  little 
explored  property,  and  that  is  the  property  that  the  Doppler  shift  is  direction 
dependent.  This  may  be  usee  as  a  means  of  direction  discrimination.  Since  the 
MTfl  is  a  noncontact  projector,  it  can  moue  at  speeds  that  would  be  quite 
impossible  for  sources  that  haue  to  be  in  physical  contact  with  the  medium. 
Therefore,  the  range  of  realizable  Doppler  shifts  can  be  extremely  large.  This 
should  lead  to  a  high  angu'ar  resolution  capability.  The  uariation  of  Doppler 
shift  with  direction  is  giuen  by  Eq.  (Ill  for  the  constant  uelocity  case.  For  the 
general  case,  it  is  necssary  to  obtain  a  numerical  solution  using  Eq.  (5). 

fin  example  is  shown  in  Fig.  5,  where  the  carrier  frequency  of  the  signal 
from  an  NTR  has  been  calculated  at  uarious  azimuth  angles.  The  MTfl  is 
tr aueling  at  M=0.9  along  the  x  axis.  Along  the  forward  direction,  the  signal 
frequency  is  20  kHz  due  to  the  Doppler  shift,  but  the  frequency  rapidly  falls  off 
away  from  the  forward  direction,  fit  24°  from  the  forward  direction,  the 
frequency  has  already  dropped  to  10  kHz.  With  a  signal  bandwidth  of  2  kHz, 
there  are  fiue  resolution  cells  in  the  24°  sector.  There  is  a  left-right  ambiguity 
which  needs  to  be  resolved,  but  that  may  be  resolued  by  using  an  asymmetrical 
path  or  by  changing  the  direction  of  movement  from  ping  to  ping. 
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Ulli.  FINDINGS 

1.  The  acoustic  signal  is  directly  a  function  of  the  modulation  of  the 
optical  power.  The  total  optical  energy  has  no  direct  effect  but  it 
must  be  sufficiently  large  to  allow  the  modulation  components  to 
exist.  This  is  generally  true  except  when  the  optical  input  is  an 
impulse,  or  if  the  source  ueiocity  component  in  the  direction  of  the 
receiver  is  equal  to  the  sound  velocity. 

2.  When  the  optical  input  is  an  impulse,  or  if  the  source  velocity  in  the 
direction  of  the  receiver  is  equal  to  the  sound  velocity,  the  acoustic 
signal  is  a  function  of  the  total  optical  energy  and  the  dimensions  of 
the  source  only.  Rny  modulation  of  the  optical  power  has  no  effect 

3.  The  moving  thermoacoustic  source  is,  at  any  speed,  neither  more  nor 
less  efficient  than  the  stationary  source  of  the  same  dimensions  and 
optical  energy. 

4.  The  impulse  train  is  the  most  efficient  optical  power  waveform  for 
the  thermoocoustic  generation  of  signals. 

5.  The  efficiency  upper  bound  is  extremely  low  for  practical  values  of 
optical  energy  and  carrier  frequency. 

6.  The  direction  dependent  Doppler  shift  of  the  acoustic  signal  gives  the 
moving  thermoacoustic  source  a  unique  property  which  cannot  be 
reproduced  by  conventional  transducers  and  which  is  potentially 
useful  for  direction  finding  without  having  to  use  a  beemformer. 

7.  As  a  laboratory  tool,  where  propagation  distances  are  short  and 
operating  frequencies  are  high,  the  Tfl  can  be  a  very  useful  tool 
particularly  because  of  Its  noncontact  property  and  wideband 
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capability. 

8.  For  sonar  applications,  a  more  efficient  optoacoustic  Conner sion  ^ 

process  is  required. 
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IK.  DISCUSSION 


These  findings  contrast  sharply  mith  some  of  the  claims  made  for  the 
efficiency  of  the  MTfl.  There  are  mainly  tioo  reasons. 

(1)  lilhen  comparing  the  efficiencies  of  transducers,  it  is  necessary  to 
define  the  desired  operating  conditions  and  then  measure  their  efficiencies 
under  the  seme  conditions.  To  make  direct  comparisons  of  efficiencies 
measured  under  different  conditions  mould  be  analogous  to  comparing  apples 
and  oranges.  One  important  aspect  of  the  operating  conditions  is  the  acoustic 
signal  itself.  For  conuentionai  transducers,  ualid  comparisons  can  be  made  if 
the  input  signals  and  the  other  releu  ant  operating  conditions  are  the  same, 
because  the  output  signal  is  normally  a  reasonable  replica  of  the  input  signal. 
This,  homelier,  is  not  the  case  for  MTfis,  mhich  can  moue  at  speeds  uery  much 
higher  and  hence  generate  far  greater  Doppler  distortions  than  conuentionai 
contact  transducers  can  ever  achseue.  Therefore  a  conscious  effort  must  be 
made  to  generate  the  some  acoustic  signal  before  any  meaningful  efficiency 
comparisons  can  be  made.  In  some  instances,8  mhiie  comparing  conuersion 
efficiencies,  the  acoustic  signal  mas  adorned  to  change  mith  the  uelocity  of  the 
MTfl  mhiie  the  optical  (input)  signal  mas  kept  constant  and  it  mas  then 
concluded  that  the  efficiency  of  an  MTfl  mas  higher  than  that  of  a  stationary 
Tfl.  This  conclusion  is  erroneous,  because  the  conuersion  efficiency  for 
generating  the  same  acoustic  signal  mas  not  compared.  The  correct  conclusion 
should  heue  been  that  the  efficiency  uoried  as  a  function  of  the  acoustic  signal 
used.  This  mistake  is  easy  to  moke  since,  os  stated  oboue,  certain 
assumptions  that  one  normady  takes  for  granted  ore  not  uabd  for  mouing 


sources. 


(2)  m  the  case  of  the  claims  by  Bozhkov,3*4  the  argument  is  reproduced 
as  follows. 

"The  efficiency ...  increases  directly  as  the  laser  intensity.  For  a  stationary 
source, ...  laser  intensity  is  limited  by  boiling  of  the  liquid  in  the  heating 
zone.  In  the  case  of  the  mouing  source  the  heat  absorbed  in  the  liquid  is 
distributed  in  space,  whereupon  the  ultimate  oalue  of  the  light  intensity 
and,  hence,  the  ultimate  efficiency  of  energy  conversion  is  greater ...  than 
for  a  stationary  source." 

This  claim  is  refuted  as  follows.  First  of  an,  it  is  the  energy  density  that  is 
limited  by  the  boiling  point  and  not  the  intensity.  However,  in  the  case  of 
water,  there  is  an  ween  hue  to  get  as  close  to  the  boiling  point  as  possible 
without  actually  causing  the  liquid  to  boil  because  the  coefficient  of  expansion, 
H,  which  directly  governs  the  conversion  efficiency,  increases  with 
temperature.  Furthermore,  should  the  peak  laser  energy  density  be  greater 
than  that  which  is  consistent  with  maximum  efficiency  without  boiling,  it  is 
possible  to  reduce  and  control  it  by  widening  the  laser  beani  ui  a  way  that  wilt 
not  increase  the  diffraction  loss  in  the  desired  signal  direction,  thereby 
obtaining  a  directional  acoustic  beam  in  the  same  process. 

Rlthough  the  MTfl  is  no  more  efficient  than  the  stationary  Tfl,  it  does 
haue  one  significant  advantage.  The  Doppler  shift  of  the  signal  is  direction 
dependent.  This  property  may  be  used  for  direction  finding.  Since  it  is  easy  to 
achieve  extremely  high  velocities,  much  higher  than  that  achievable  with 
transducers  that  have  to  be  in  contact  with  the  water,  the  potential  angular 
resolution  may  be  quite  impressive. 


H.  CRffiCRL  ISSUES 


fl.  Efficiency  Enhancement 

Clear  kj,  the  efficiency  of  the  thermoacoustic  process  is  too  low  to  be 
useful  in  sonar  applications  except  at  very  short  ranges.  R  radically  different 
process  that  has  the  advantages  of  remote  deployment  and  noncontact  with 
the  ivater  is  needed.  The  so-caned  nonlinear  optoacoustic  process9  is  ivorth 
investigating.  High  efficiencies  (up  to  10%)  have  been  claimed  for  it,  although 
this  may  be  somewhat  over  optimistic.  It  has  been  investigated  by  several 
researchers,9"13  but  the  resists  published  so  far  are  not  applicable  to  the 
sonar  application. 

Rs  a  first  attempt,  a  simple  blast  model  should  be  used  to  predict  the 
acoustic  output  of  the  optoacoustic  process.  R  possible  candidate  blast  model  is 
that  of  R.  N.  Pirri1^  who  successfully  used  it  to  model  the  momentum  transfer 
between  a  high  power  laser  blast  on  the  surface  of  a  solid.15  R  more 
compbcated  model  has  been  reported  but  it  is  essentially  a  one-dimensional 
me  iel 16  deueloped  for  solid  surfaces. 

B.  Planned  Work 

The  blast  model  needs  to  be  developed  and  verified  experimentally 
unC’!, r  controlled  laboratory  conditions.  Reliable  riuontitotive  measurements  of 
the  optoacoustic  sound  generation  process  have  so  far  not  been  found  in  the 
published  literature.  The  reported  works  have  been  either  qualitative10,11  or 


concerned  with  Impulse  generation,1^’13  inhere  most  of  the  energy  is  spread 
over  a  band  of  seueral  megahertz,  with  one  exception,17  inhere  the 
experimental  conditions  appeared  to  be  imprecise.  fi  Q-sinitch  driner  mill  be 
purchased  in  the  fail  of  1984  to  enable  the  laser  to  deliuer  short  pulses  of 
sufficient  intensity  to  generate  blasts.  Experiments  mill  be  carried  out  to 
measure  the  acoustic  output  and  its  spectrum.  From  the  experimental  results, 
a  more  realistic  model  of  the  optoacoustic  connersion  process  may  be  obtained. 
Then  a  reliable  estimator  of  system  performance  nersus  laser  power 
requirements  can  be  constructed  for  further  evaluation. 


RPPEMMHfl 

BflSIC  THERMOftCOUSTIC  GENERATION  THEORY 


The  acoustic  pressure  may  be  deriued  from  the  inhomogeneous  waue 

equation 

Apa  -  (a2pa/3tV2  = -3q(t-r/c)/3t .  (A.1) 

For  a  thermoacoustic  source,  the  source  function  q(t)  is  giuen  by 

q(t)  =  -  V.l(x,y,z,t)  K/Cp .  (fl.2) 

For  a  stationary  TA  whose  axis  passes  through  the  origin  slanted  at  an  angle  81 
reiatiue  to  the  uertical,  directed  at  an  azimuth  angle  0^  relatiue  to  the  x  axis, 
and  an  intensity  profile  !$(x,y,t)  at  the  surface, 


V.l(x,y,z,t)  *  A  Uuj.t^ji  exp(-|*2/cos(0 ,  (fl.3) 

where 

Xy,  *  x  -  z  costO^  sln(81) 
and 

«  y  -  z  coshty  sintB^ . 

It  is  assumed  that  the  medium  reacts  to  the  thermal  input  instantaneously. 
Taking  the  Fourier  transform  of  Eq.  (A.1), 

APaM  ♦  k2  Pa(w)  -  AjicH/Cp  expf-jiz/costBj))  k  ls(  y^w) ,  (R.4) 

49 


J 


where 

k  =  w/c , 

Pa(<*>)  =  pa(t)  expOwt)  dt , 

and 

ls(  Xy,,  y^w)  =  _<*/“  l^x^  y^t)  exp(iwt)  dt . 

Solukig  Eq.  (fl.4)  for  the  farfield  solution,  the  following  expression  is  obtained 
for  the  acoustic  pressure  Fourier  transform, 

PaM  -  BjJtcH/Cp  -oJt0.cJC°-cJCO  -G(r,,r0)  expf-jiz/cosfS^f 

k  l$(  xw,  yw,w)  dxdydz,  (fl.4) 

where 

^l'^,)  *  exp(-ilcj  -  r  1  -  ar)/(4n  |  -  cD 

*  Pau^  expHki  f^,  -  r1  -  ar)/(4ir  |  tg  -  rl) , 

C"  is  the  mirror  image  of  c'  in  the  surface,  and 

Paw  ls  ttie  acoustic  reflection  coefficient  of  the  water-air  surface;  it 

is  assumed  to  be  equal  to  -1, 
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The  geometry  is  shown  in  Fig.  7.  Rewriting  Eg.  (fl.4), 


Pgtw)  =  k  Ca  Da  U)0M  exp(  ikr0  -  ar0  )  /r0 


(R.5) 


where 

U)Qm  =  _m!m  w0(t)  exptiwt)  dt 

=  r-mr  !s(  y^w)  dx  dy , 

Ca  =  RKc/(4TfCp) , 

and 

08  *  -oJ^-cJ^-cJ00  (|»  /cus  8t)  exp(-j*z/cos  8,) 

(  -GCCxq)  )expHkr0+  oxQ )  ls(  x^  y^w)/  uyw)  dx  dy  dz 


For  a  uertical  thermoacoustic  array,  the  ualue  of  8^  is  zero  and  the 
expression  for  0a  reduces  to 

Da  -  ,Mr >1  espl-uz) 

i  -Gir'^g)  )expHkr0*  ar0 )  ls(  x^  y^wy  UJ0M  dx  dy  dz . 


FIGURE  7 

THE  GEOMETRY  OF  A  GENERAL  STATIONARY  THERMOACOUSTIC  ARRAY 
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APPENDS  B 

MfiHMUM  OPTICA  POWER  FOURIER  TRftNSFORM 
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Giuen  a  real  waveform,  w(t),  and  its  Fourier  transform,  LiMu>),  and 
given  that  wit)  is  greater  than  or  equal  to  zero,  the  upper  hound  of  the  ratio 
ID10)  is  determined  as  follows. 

Since  w(t)  >  0,  then  w(t)  may  be  expressed  as  the  square  of  another 
function  z(t), 

wW-fzWf2,  (B.1) 

where  ztt)  is  a  complex  wane  form. 

Eq.  (B.1)  may  be  rewritten  as 

w(t) 3  z(t)  z*(t) .  (B 2) 

Let  Z(o>)  be  the  Fourier  transform  of  ztt).  From  Eq.  (B^),  it  follows  that  llMw)  is 
t.ie  convolution  of  Ztoi)  and  2(-0)L  Since  the  convolution  of  any  function  with 
the  reverse  of  itself  ts  equivalent  to  autocorrelation,  IDtw)  is  therefore  the 
autocorrelation  function  of  2(0)).  Thus 

UXw)  -  2(0),)  Z*(<a>,-<->)  do), .  (B.3) 

For  any  bounded  function,  the  peak  of  the  autocorrelation  function  is  real  and 
positive  and  always  occurs  at  zero  dwptacement.  Therefore, 

(WO)  >  I  U)M .  (B.4) 
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